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We  study  complex  processes  whose  evolution  in  time  rests  on  the  occurrence  of  a  large  and  random  number  of  events.  The  mean 
time  interval  between  two  consecutive  critical  events  is  infinite,  thereby  violating  the  ergodic  condition  and  activating  at  the  same 
time  a  stochastic  central  limit  theorem  that  supports  the  hypothesis  that  the  Mittag-Leffler  function  is  a  universal  property  of  nature. 
The  time  evolution  of  these  complex  systems  is  properly  generated  by  means  of  fractional  differential  equations,  thus  leading  to  the 
interpretation  of  fractional  trajectories  as  the  average  over  many  random  trajectories  each  of  which  satisfies  the  stochastic  central 
limit  theorem  and  the  condition  for  the  Mittag-Leffler  universality. 


1.  Introduction 

The  fractional  calculus  has  developed  in  a  number  of  signif¬ 
icant  ways  in  the  recent  past.  Sokolov  et  al.  [1]  maintain  that 
this  calculus  was  restricted  to  the  field  of  mathematics  until 
the  last  decade  of  the  twentieth  century,  when  it  became  very 
popular  among  physicists  as  a  powerful  way  to  describe  the 
dynamics  of  a  variety  of  complex  physical  phenomena.  For 
example,  anomalous  diffusion  was  described  using  fractional 
diffusion  equations  [2,  3];  viscoelastic  materials  were  mod¬ 
eled  using  fractional  Langevin  equations  [4];  and  complex 
dynamic  systems  could  be  governed  using  fractional  control 
[5].  In  the  last  decade  the  concept  of  fractional  dynamics 
has  gained  further  attention  in  the  statistical  and  chemical 
physics  communities  [6].  Fractional  differential  equations 
have  also  been  successfully  applied  to  neural  dynamics  [7,  8] 
and  ecology  [9]  as  well  as  to  traditional  fields  of  engineering 
[10, 11]  namely  [12, 13]. 

Of  particular  interest  to  the  authors  is  the  growing  liter¬ 
ature  on  extending  systems  of  nonlinear  dynamic  equations 
having  strange  attractor  solutions  to  fractional  nonlinear 
equations.  Such  extensions  were  typically  made  by  replac¬ 
ing  integer-valued  derivatives  by  fractional  derivatives;  for 
example,  in  the  Lorenz  system  [14-16],  in  the  chaotic  rigid 


body  motion  of  gyros  [17],  in  Hopfield-type  neural  networks 
[8],  and  in  the  immune  model  of  HIV  infection  [18],  to 
name  a  few.  These  replacements  were  made  in  attempts  to 
incorporate  dynamic  mechanisms  thought  to  be  important 
that  could  not  be  captured  by  the  traditional  models,  for 
example,  complexity  in  the  form  of  memory  in  time  and  non¬ 
locality  in  space.  The  results  of  extending  these  nonlinear 
models  has  been  to  apparently  introduce  dissipation  into 
the  dynamics  such  that  the  solution  on  the  strange  attractor 
collapses  to  that  of  a  stable  fixed  point.  The  appearance 
of  a  fixed  point  is  interpreted  to  be  a  consequence  of  an 
induced  dissipation  mimicking  the  complexity  modeled  by 
the  fractional  derivatives. 

Herein  we  provide  an  alternative  interpretation  of  these 
extensions  that  involves  the  notion  of  a  fluctuating  trajectory 
and  interpreting  the  fractional  models  as  averages  over  an 
ensemble  of  these  trajectories  on  the  strange  attractor.  This 
view  is  consistent  with  one  proposed  as  an  extension  of  con¬ 
servative  Hamiltonian  systems  to  fractional  systems  [19] .  This 
generalization  of  classical  mechanics  is  based  on  a  random¬ 
ization  of  chronological  or  clock  time  in  the  traditional  phase 
space  using  the  notion  of  operational  time  and  subordination. 
Without  going  into  the  details  of  the  extension  of  Hamilton’s 
equations  of  motion  to  fractional  form  it  suffices  to  note  that 
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fractional  derivatives  in  chronological  time  are  interpreted 
as  averages  of  a  particle’s  displacement  and  momentum  over 
the  fluctuating  operational  time.  Consequently,  the  dynamics 
of  a  single  fractional  harmonic  oscillator,  for  example,  is 
considered  to  be  an  average  over  an  ensemble  of  harmonic 
oscillators.  Stanislavsky  [19]  emphasizes  that  each  oscillator 
differs  slightly  from  every  other  oscillator  in  frequency 
because  of  subordination.  The  phase  space  trajectories  rather 
than  being  level  energy  curves  instead  spiral  into  the  origin 
and  the  fractional  oscillator  “demonstrates  a  dissipative  pro¬ 
cess  stochastic  by  nature.” 

For  the  more  general  dynamical  systems  considered 
herein  there  is  no  Hamiltonian  with  which  to  generate  the 
equations  of  motion.  Recall  that  a  strange  attractor  requires 
the  system  to  dissipate  energy.  Therefore  the  approach  taken 
here  is  stochastic  rather  than  dynamic  and  requires  the  devel¬ 
opment  of  a  new  form  of  the  central  limit  theorem  (CLT), 
one  that  is  compatible  with  the  fractional  calculus.  Section  2 
connects  the  familiar  Mittag-Teffler  function  (MLF)  solution 
to  a  Caputo  fractional  differential  equation  [20]  to  what  we 
have  named  the  stochastic  central  limit  theorem  (SCLT). 

At  the  most  elementary  level  the  familiar  relaxation  rate 
equation  is  replaced  by  its  fractional  form 


(t) 

dta 


-A“'F  (f) 


(1) 


and  here  we  interpret  the  fractional  derivative  to  be  of  the 
Caputo  form 


s  i  r  dt'  w  (fl) 

dta  I  (a  -  ri)  Jo  (t  -  t')a+1~n  V  ’ 


(2) 


where  the  superscript  is  the  derivative  of  integer  order  n  such 
that  n- 1  <  a  <  n  and  since  a  <  1  we  have  n  =  1.  The  solution 
to  (1)  is  the  MLF  [21, 22] 


^ML  (t)  =  Ea  (-(Af)“) 
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(4) 

as  t  — >  0, 

ast  — »  oo. 

(5) 


Metzler  and  Klaffer  [23]  exploited  this  connection  between 
the  stretched  exponential  and  inverse  power-law  (IPL)  distri¬ 
butions  by  interpreting  the  MLF  as  a  survival  probability  and 
thereby  establishing  a  bridge  between  the  advocates  of  these 
two  distinct  forms  of  the  survival  probabilities  as  important 
signs  of  complexity. 

In  the  Materials  and  Methods  section  we  show  that 
TML(f)  is  universal,  as  had  been  advocated  by  Gorenflo  and 
Mainardi  [24]  and  others  [25]  in  the  same  sense  as  the 
limit  distributions  in  the  CLTs  of  Laplace  (CLT)  and  Levy 
(generalized  CLT)  but  here  the  universality  is  shown  to  be 


a  consequence  of  the  SCLT.  In  Section  2.1  the  form  of  the 
SCLT  is  set  up  and  in  Section  2.2  the  theorem  is  proven 
using  a  scaling  argument.  As  a  consequence  of  the  SCLT  we 
show  in  Section  2.3  by  introducing  fluctuating  trajectories, 
interpreted  as  a  proper  representation  of  complex  processes 
with  memory,  that  the  MLF  is  universal.  The  universality 
is  a  consequence  of  a  subordination  process.  In  the  Results 
and  Discussion  section  the  connection  between  the  Caputo 
fractional  derivative  and  TML(f)  together  with  the  fractional 
trajectories  of  Section  2,  through  subordination  result  in  the 
fractional  derivative  being  interpreted  as  an  average  over 
infinitely  many  random  trajectories.  In  Section  4  we  draw 
some  conclusions. 


2.  Materials  and  Methods 

In  the  late  1980s  there  was  remarkable  activity  in  the  develop¬ 
ment  of  the  theory  of  random  summations,  namely,  the  case 
where  the  number  of  summands  is  itself  a  random  variable 
[26].  We  adopt  the  term  stochastic  limit  theorem  rather  than 
random  summation  to  emphasize  that  we  depart  from  the 
exemplary  Poisson  condition  of  [26] .  The  MLF  is  generated 
by  a  fluctuating  number  of  events  m  with  fluctuations  as 
large  as  the  mean  value  (m).  The  main  difference  between 
the  traditional  and  the  SCLT  is  that  the  former  rests  on  the 
sum  of  a  fixed  number  m  of  fluctuations,  and  on  the  rescaling 
procedure  to  use  for  m  — >  oo.  The  SCLT  is  based  on 
keeping  fixed  the  probability  Ps  of  detecting  an  event,  that 
is,  the  probability  that  an  event  is  visible  in  an  experiment. 
Each  value  of  Ps  generates  a  sequence  of  m  elementary 
laminar  regions  but  only  one  visible  event  at  the  end  of 
the  last  laminar  region.  The  SCLT  focuses  on  the  interval 
between  two  consecutive  visible  events  and  adopts  a  rescaling 
procedure  to  compensate  for  the  incomplete-measurement- 
induced  survival  probability  enhancement  [27].  We  show  that 
all  waiting-time  pdf’s  generating  a  non-integrable  survival 
probability  as  a  consequence  of  the  SCLT  yield  i//ML(t).  This 
proof  leads  us  to  conclude  that  vI/ML(t)  is  universal. 


2.1.  Stochastic  Central  Limit  Theorem.  The  concept  of  survival 
probability  is  connected  to  the  stochastic  perspective  of  a 
complex  system  generating  events  in  time.  The  time  interval 
between  two  consecutive  events  (laminar  region)  is  assigned 
the  values  ±1,  according  to  a  coin  tossing  prescription  [28]. 
At  time  t  =  0  the  system  is  prepared  by  selecting  all  the 
realizations  with  an  event  occurring  at  that  time,  with  ensuing 
positive  laminar  regions.  As  a  consequence  the  probability 
that  no  event  occurs  up  to  time  t,  denoted  as  T(f),  is  properly 
termed  a  survival  probability,  and  the  function 


f(t)  =  - 


dW  ( t ) 
dt 


(6) 


is  the  waiting-time  probability  density  function  (pdf).  We 
adopt  the  symbols  xFML(t)  and  i//ML(t)  to  denote  the  MLF 
survival  probability  and  the  corresponding  waiting-time  pdf, 
respectively.  Physical  examples  of  y/ML(t)  generated  by  the 
cooperative  interaction  of  many  units  can  be  found  [27,  29], 
with  the  important  observation  that  the  stretched  exponential 
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regime  becomes  more  extended  if  the  probability  of  gener¬ 
ating  a  visible  cooperative  event  decreases  as  discussed  in 
Section  2.2. 

We  assume  that  the  time  interval  between  two  consecu¬ 
tive  critical  events  generated  by  the  complex  system  under 
study  is  given  by  the  waiting-time  pdf 

i//(r)  oc  with  1  <  p  <  2.  (7) 

The  corresponding  cumulative  distribution  'F(r)  has  the  form 
lim  T  (t)  oc  — ,  with  a  =  u  -  1  <  1.  (8) 

r^oo  w  T<*  r 


thereby  yielding  for  the  subsidiary  function  in  (10) 


E(u)  = 


u  \  1 

T0  /  1  +  (u/A0) 


(15) 


satisfying  the  condition  of  (11).  In  other  words,  the  properties 
of  (10)  and  (11)  are  fulfilled  by  all  the  waiting-time  pdf’s  with 
the  scale-free  condition  of  (7) .  We  now  show  that  the  waiting¬ 
time  pdf  corresponding  to  the  sum  of  a  large  numbers  of 
times  each  of  which  is  generated  by  the  generic  pdf  y/(t)  of 
(7),  not  necessarily  of  the  MLF  type,  are  MLF  waiting-time 
pdf’s. 


The  origin  of  this  condition,  usually  interpreted  as  a  mani¬ 
festation  of  complexity,  can  either  be  the  anomalous  nature 
of  the  dynamics  under  investigation  [30,  31]  or  the  condition 
of  criticality  [32].  In  the  former  case  the  property  described 
by  (7)  can,  for  example,  be  the  consequence  of  diffusing 
molecules  being  trapped  for  long  times  in  wells  with  a 
random  distribution  of  depths.  In  the  latter  and  less  well 
known  situation  the  emergence  of  temporal  complexity  is 
due  to  the  cooperative  action  of  many  interacting  units.  At 
the  onset  of  the  cooperation-induced  phase  transition  from 
disorder  to  order,  the  mean  field  fluctuates  and  its  non¬ 
stationary  waiting-time  pdf  corresponds  to  an  IPL  y/(t)  [33]. 

We  adopt  for  the  Laplace  transform  of  the  time  function 
/(f)  the  following  notation: 


2.2.  Imperfect  Detection  of  Events.  To  make  this  demonstra¬ 
tion  as  clear  as  possible  and  at  the  same  time  provide  an 
intuitive  understanding  of  the  SCLT,  let  us  imagine  that  the 
detector  used  to  monitor  the  events  produced  by  (7)  is  not 
very  accurate  and  that  the  probability  of  perceiving  these 
events  is 

PS<  1-  (16) 

As  a  consequence  of  this  imperfection  a  time  t  between  two 
consecutive  visible  events  is  the  sum  of  m  elementary  times 
derived  from  the  condition 

P  (m)  =  Ps(l  -  Ps)m,  (17) 


rOO 

f  (u)  =  £{f  (t);u}  =  dt  exp  (-ut)  f  (f) .  (9) 

Jo 

It  is  important  to  stress  that  to  satisfy  the  long-time  limit  of 
(7)  the  Laplace  transform  of  i//(r)  has  the  functional  form 


which  is  the  probability  that  the  first  m  events  after  the  initial 
preparation  event  are  not  visible  while  the  (m  +  l)th  event  is 
visible.  For  Ps  — >  0, 

P  (m)  =  Ps  exp  ( -mPs ) ,  (18) 


/  (u)  =  1  - 


+  H  (u) , 


(10) 


thereby  implying  that  the  standard  deviation  is  on  the  same 
order  as  the  mean 


with  the  condition  on  the  subsidiary  function  H (u) 


lim  A 

u  — >  0 


a 

0 


=  o. 


(ii) 


Therefore  the  subsidiary  function  must  vanish  more  rapidly 
than  ua  as  u  — >  0. 

Note  that  the  Laplace  transform  of  the  MLF  survival 
probability  given  by  (3)  is  [20, 22] 


^ml  (M)  - 


ua  +  A“ 


(12) 


and  the  relation  to  the  Laplace  transform  of  the  waiting-time 
pdf  is 


^ml  (M)  ~ 


1  ~  V^ml  («) 
u 


(13) 


so  that  with  a  little  algebra  we  obtain 


((m2)  -  <m>2) 

(m>2 

Thus,  the  condition 


(19) 


,  ,  1 

(m)  =  - - »  oo  (20) 

Ts 

of  the  SCLT  is  quite  different  from  the  condition  m  — »  oo  of 
the  Gauss  and  Levy  CLTs,  since  m  has  very  large  fluctuations 
around  (m)  in  the  traditional  argument.  To  better  understand 
the  new  theorem,  we  note  that  the  probability  of  generating 
at  time  t  an  event  that  is  the  last  of  a  sequence  of  m 
events  occurring  at  earlier  times,  i//m(f),  does  not  satisfy  the 
condition  of  generating  a  MLF  stable  form  for  m  — »  oo. 
Llowever,  the  waiting-time  pdf  of  the  time  intervals  between 
visible  events,  does,  for  Ps  — »  0.  The  pdf  of  finding  a  visible 
event  a  time  interval  t  after  an  earlier  visible  event  is  given  by 


V^ML  (M)  - 


1 

(w/A0)“  +  1  ’ 


W  (1)  =  I  PS(1  -  PS)mfm+,  «  • 

m=0 


(14) 


(21) 
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Implementing  the  well  known  property  for  renewal  processes  It  is  straightforward  to  show  that  inserting  fv(u)  of  (23)  into 

[22,  28]  (30)  yields  a  result  equivalent  to 


fm(u)  =  ty(u)]m  (22) 

we  obtain  from  the  Laplace  transform  of  (21)  after  a  little 
algebra 


VV  (u)  = 


_ yM _ 

1-((1  -Ps)/Ps)  (£(«)-!)' 


(23) 


Using  the  partitioning  of  (10)  we  rewrite  (23)  in  the  more 
convenient  form 


-  ,  .  _  _ 1  -  (m/A0)“  +  5  (u) _ 

fv  U>  “  1  +  ((1  -  Ps)  /Ps)  (u/ A0)“  -  ((1  -  Ps)  /Ps)  H  ( u ) ' 

(24) 

To  fulfill  the  limiting  condition  of  (11)  the  slowest  contribu¬ 
tion  to  E(u)  must  be 


^slowest  («)  =  kua+e,  (25) 

withe  >  0.  Rescaling  the  Laplace  variable  u  with  the  detection 
probability 


u 


=  uPl'a, 


(26) 


transforms  (24)  into 


Vv(w')  = 


that  for  Pq 


l-Ps(u'/X0)a  +  E(u'PsUa) 


1  +  (1  -  Ps)  («'M 0)“  -  ((1  -  ^s)  /i’s) 5  (w'Ps1/a) 

(27) 

— »  0  reduces  to 


VV  («')  =  “ — 77 

1  +  (l/ 


(u'/A0)a 


(28) 


In  fact,  the  rescaled  slowest  contribution  to  E(u )  is  pro¬ 
portional  to  Ps1+e^“  thereby  making  the  contributions  of 
H  (u)  vanish  for  Ps  — >  0  in  both  the  numerator  and  the 
denominator  of  (27). 

We  see  that  after  rescaling  frv(u)  coincides  with  (14)  and 
consequently  its  inverse  Laplace  transform  is 


VV  (0  = 


dEa(-(X0t)a ) 


dt 


(29) 


Consequently,  the  survival  probability  of  the  visible  events  is 
the  MLF.  This  is  the  essence  of  the  SCLT. 


2.3.  Fluctuating  Trajectories.  To  use  the  SCLT  to  interpret 
the  fractional  trajectories  as  averages  over  infinitely  many 
stochastic  realizations,  it  is  convenient  to  get  vPML(t)  from 
the  survival  probability  of  the  time  interval  between  two 
consecutive  visible  events.  We  note  from  (13)  that 


TV  (u)  = 


1  -  VV  (u) 

u 


(30) 


OO  ,t 

(t)  =  I|  (0  v(t-  0  (1  -  PsY>  (31) 

n= 0  Jo 


which  has  the  well  known  Montroll-Weiss  continuous  time 
random  walk  (CTRW)  structure  [34] .  The  Laplace  transform 
of  TV(0  from  (31)  reads 


WV(U)  = 


1 

u  +  O  (u)  Ps 


(32) 


where  &(u)  is  the  Laplace  transform  of  the  Montroll-Weiss 
memory  kernel  defined  by 


<J>  ( u )  = 


u\j/  ( u ) 

1  -  \j/(u) 


(33) 


Note  that  by  using  inverse  Laplace  transforms  it  is  straight¬ 
forward  to  establish  that  (31)  is  equivalent  to 


jVy  (; t )  =  -Ps  J‘  dt’O  (t’)  Wv  ( t  -  t’)  ,  (34) 


whose  time  convolution  structure  justifies  the  adoption  of  the 
term  memory  kernel  for  O(t). 

Due  to  the  equivalence  between  (21)  and  (31)  we  can  use 
the  Ps  scaling  argument  to  immediately  conclude  that 


1 

u'  +un~a  A“’ 


(35) 


thereby  demonstrating  why  the  MLF  is  ubiquitous  in  data  sets 
since  it  is  the  distribution  of  visible,  that  is,  measurable  events. 

To  shed  further  light  into  the  SCLT,  let  us  notice  that  the 
condition  Ps  — »  0  has  the  effect  of  turning  n  into  a  virtually 
continuous  time  r,  (1  -  Ps)n  into  exp(-Psr)  and  T rv(t)  of  (31) 
into  Ea(-Ps(X0t)a).  Thus  the  MLF  survival  probability  with 
A  =  (Ps)1/aA0  is  the  counterpart  in  real  time  t  of  the  ordinary 
exponential  function  exp(-Psr), 

rOO 

Ea  (-Fs(vr)  =  Jo  drp^  (t, r)  exp  (-Pst)  (36) 


and  p®(t,  r)  is  the  pdf  of  times  t  corresponding  to  the 
continuous  time  r.  A  straightforward  way  of  deriving  (36) 
from  (31)  rests  on  observing  that  the  Laplace  transform  of 
(31),  namely,  (32),  can  be  interpreted  as  a  double  Laplace 
transform  p ^  ( u ,  s)  with  s  =  Ps.  By  inverse  Laplace  trans¬ 
forming  p<s)  (u,  s  =  Ps)  with  respect  to  u  we  obtain  (36)  in 
accordance  with  [35-37].  The  time  t  can  be  interpreted  as  a 
diffusing  position  of  a  random  walker  that  keeps  jumping  in 
the  same  direction.  It  is,  in  fact,  an  asymmetric  Levy  process 
[35-37].  We  are  therefore  led  to  interpret  the  SCLT  as  a 
consequence  of  the  generalized  CLT  of  Levy.  The  condition 
n  — »  oo  generates  the  Levy  stable  pdf  and  the  sum  over 
infinitely  many  Levy  processes  weighted  by  the  exponential 
function  exp(-Psr)  generates  the  MLF  survival  probability, 
in  accordance  with  earlier  results  [38,  39]. 
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3.  Results  and  Discussion 

The  SCLT  is  the  first  important  result  of  this  paper.  The  tra¬ 
ditional  CLT  yields  the  Normal  distribution  as  the  limit  pdf. 
The  generalized  CLT  developed  by  Levy  yields  the  a-stable 
distribution  as  the  limit  pdf.  Finally  the  SCLT  presented 
herein  yields  the  MLF  as  the  limit  pdf.  On  the  basis  of  the 
SCLT  we  established  a  physical  interpretation  of  fractional 
trajectories  that  had  not  been  previously  considered  in  the 
literature  for  non-conservative  dynamic  systems.  The  MLF 
survival  probability  can  be  interpreted  as  the  realization  in 
the  continuous  time  t  of  the  traditional  exponential  relaxation 
process 


dT  (r) 

dr 


-Wt) 


(37) 


in  operational  time.  The  probabilistic  structure  of  (36)  rep¬ 
resents  the  sum  over  infinitely  many  such  random  relaxation 
processes,  which  turns  out  to  be  equivalent  to  the  MLF  sur¬ 
vival  probability.  In  the  same  way  a  fractional  trajectory,  not 
necessarily  dissipative,  was  proven  to  be  a  sum  over  infinitely 
many  stochastic  trajectories  that  for  historical  reasons  we  call 
Montroll-Weiss  trajectories. 


The  numerical  procedure  replaces  the  continuous  one  with 
a  discrete  time  as  would  be  appropriate  in  a  numerical 
algorithm.  Thus  we  can  write  for  the  nth  iterate  of  this 
equation,  using  the  notation  V(n)  =  V(nAr): 

V  (n)  =  (l  -  AtF)"v  (0) .  (41) 

Note  that  the  numerical  solution  of  (38)  realized  by  means 
of  the  prescription  (41)  is  an  extension  to  a  generic  trajectory 
of  the  exponential  relaxation  (1  -  Ps)n.  For  this  reason  the 
condition  of  perfect  integration  Ar  — »  0  corresponds  to 
Ps  —>  0,  thereby  establishing  a  connection  with  the  SCLT. 

A  stochastic  Montroll-Weiss  trajectory  is  obtained  by 
assuming  the  transition  V(«)  — >  V(»  +  1)  is  a  crucial  event 
occurring  at  the  time  t(n)  =  +  •  •  •  rn,  where  the  r;  are 

random  times  with  the  waiting-time  pdf  of  (7),  not 
necessarily  identical  to  t//ML(r).  We  assign  to  V(f)  the  value 
V(n)  with  n  fulfilling  the  condition  t(n)  <t<  t(n  +  1). 

It  is  important  to  restate  the  difference  between  the  Levy 
CLT  and  the  SCLT.  As  shown  in  [30],  the  pdf  p(s\t,r) 
is  an  asymmetric  Levy  process,  which  ought  not  to  be 
confused  with  the  stable  i//ML(t)  generated  by  the  procedure 
of  this  paper.  An  average  over  infinitely  many  Montroll-Weiss 
trajectories  yields 


3.1.  Subordination  and  Fractional  Dynamics.  To  relate  the 
physical  interpretation  established  in  the  previous  section  to 
the  replacement  fractional  differential  equations  we  intro¬ 
duce  the  nomenclature  of  a  fractional  trajectory.  Consider  the 
differential  equation  in  operational  time 

-y-V  (r)  =  -TV  (r) ,  (38) 

dr 

where  V  is  a  generic  multidimensional  vector  and  F  is  a 
generic  operator  either  linear  or  nonlinear  acting  on  the 
components  of  the  vector  V.  With  this  concise  notation 
we  may  describe,  for  instance,  the  Lorenz  system  [14-16], 
which  is  a  remarkable  application  of  the  fractional  calculus 
to  chaos,  with  V  =  {X,Y,Z}.  All  the  studies  using  the 
fractional  calculus  can  be  represented  by  the  notation  of 
(38).  Consequently,  in  the  nonlinear  systems  whose  dynamics 
unfold  on  a  strange  attractor  we  follow  the  tradition  and 
replace  the  integer  order  derivative  d/dr  with  the  Caputo 
fractional  derivative  of  index  a,  with  a  <  1,  to  obtain 

£v(f)  =  -rv(f).  (39) 

We  refer  to  the  solution  of  (39)  as  a  fractional  trajectory,  that 
being  the  time  trace  of  the  solution  in  phase  space  for  the 
system.  Here  we  anticipate  the  second  result  of  this  paper 
that  being  the  time  evolution  of  V (t) ,  the  fractional  trajectory, 
is  an  average  over  infinitely  many  stochastic  Montroll-Weiss 
trajectories. 

To  properly  define  a  Montroll-Weiss  trajectory,  let  us  go 
back  to  (38)  and  adopt  a  numerical  procedure  to  integrate  it 
over  the  operational  time  and  obtain 

V(t  +  At)  =  (i-AtF)V(t).  (40) 


V(t)  = 


oo  rt 

X  dt'Vn  (t')  xp(t-t')(l-  Arr)”v  (0) . 

n= 0  Jo 


(42) 


The  left  hand  side  of  (42)  is  a  Gibbs  ensemble  average  over 
infinitely  many  fluctuating  trajectories  and  should  be  denoted 
by  the  symbol  (V(f)).  However  to  stress  its  connection  with 
the  fractional  derivative  formalism  we  continue  to  use  the 
symbol  V(t).  The  Laplace  transform  of  V(f)  is 


V  (u)  =  - ^ - =V(0). 

u  +  O  (u)  Arr 

On  the  basis  of  the  SCLT,  we  conclude  that  for  At 
corresponding  to  Ps  — >  0, 


(43) 

0, 


V  (u)  = 


- =V(0). 

u  +  u1  “AJjArr 


(44) 


This  allows  us  to  replace  the  average  given  by  (42)  with  its 
fractional  differential  equivalent 

^V(f)  =  -ArA“fV(f).  (45) 

In  fact,  using  the  Caputo  fractional  derivatives  we  have 


£  - 


d“ 

df“ 


V  (f) ;  u 


=  m“V  (m)  -  ua~l  V  (0) . 


(46) 


Laplace  transforming  (45)  with  this  rule  has  the  effect  of 
yielding  (44),  thereby  establishing  that  the  fractional  differen¬ 
tial  equation  given  by  (39)  is  equivalent  to  the  mean  Montroll- 
Weiss  trajectory  of  (42),  under  the  condition  A^At  =  1. 
Note  that  At  must  be  small  enough  as  to  ensure  a  convergent 
solution  to  (38)  with  the  possible  effect  of  making  A0  so  large 
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as  to  virtually  cancel  the  stretched  exponential  regime  of  the 
MLF.  With  AqAt  =  1  the  effective  rate  A  is  determined  by 
the  eigenvalues  of  f,  which  determine  the  density  of  events, 
thereby  making  the  stretched  exponential  become  ostensible 
with  very  low  event  densities,  in  accordance  again  with  the 
property  of  incomplete  measurement  [27]. 

4.  Conclusions 

We  have  established  the  universality  of  the  MLF  survival 
probability  using  a  SCLT.  The  term  stochastic  in  SCLT 
emphasizes  the  fact  that  the  large  number  of  elementary 
laminar  times  whose  sum  generates  the  time  interval  between 
two  consecutive  visible  events  is  a  widely  fluctuating  number 
with  a  fluctuation  intensity  as  large  as  the  mean  number  of 
events  involved  in  the  process.  The  SCLT  leads  to  a  new 
perspective  of  fractional  trajectories,  which  yields  a  new 
physical  interpretation  of  their  claimed  stronger  stability 
on  the  part  of  those  that  have  replaced  the  integer  with 
fractional  derivatives.  The  dissipative  nature  of  the  fractional 
trajectories  has  to  be  interpreted  as  a  form  of  phase  decor¬ 
relation  process  rather  than  one  with  friction.  The  fractional 
version  of  the  popular  Lotka-Volterra  ecological  dynamics, 
for  instance,  generates  a  fractional  trajectory  that  is  the 
sum  over  infinitely  many  Montroll-Weiss  trajectories.  Each 
Montroll -Weiss  trajectory  is  an  ordinary  Lotka-Volterra  cycle 
in  the  operational  time  n.  Transitioning  from  the  operational 
time  n  to  the  chronological  time  t  spreads  these  trajectories 
over  the  entire  Lotka-Volterra  cycle  thereby  generating  the 
mistaken  impression  that  the  resulting  average  trajectory 
reaches  equilibrium  through  a  dissipative  process. 
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